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Abstract 

It has been known for a long time that for birth-and-death processes started in zero the 
first passage time of a given level is distributed as a sum of independent exponentially 
distributed random variables, the parameters of which are the negatives of the eigenvalues 
of the stopped process. Recently, Diaconis and Miclo have given a probabilistic proof 
of this fact by constructing a coupling between a general birth-and-death process and a 
process whose birth rates are the negatives of the eigenvalues, ordered from high to low, 
and whose death rates are zero, in such a way that the latter process is always ahead of 
the former, and both arrive at the same time at the given level. In this note, we extend 
their methods by constructing a third process, whose birth rates are the negatives of the 
eigenvalues ordered from low to high and whose death rates are zero, which always lags 
behind the original process and also arrives at the same time. 
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1 Introduction 

1.1 First passage times of birth-and-death processes 

Let X = {Xt)t>o be the continuous-time Markov process in N = {0,1,...}, started from 
Xq = 0, that jumps from x — 1 io x with hirth rate bx > ^ and from x to x — 1 with death rate 
4 > (x > 1). Let 

TN := inf{t >Q:Xt = N} {N >l) (1.1) 

denote the first passage time of N . The following result has been known at least since |KM591 
Prop. 1]. 

Proposition 1.1 (Law of first passage times) The first passage time is distributed as 
a sum of independent exponentially distributed random variables whose parameters Ai < • • • < 
A AT are the negatives of the nonzero eigenvalues of the generator of the process stopped in N. 



Older proofs of this fact are based on a calculation of the Laplace transform of by 
purely algebraic methods, see |DM09| for a historical overview. In the latter paper, Diaconis 
and Miclo gave for the first time a probabilistic proof of Proposition 11.11 by coupling the 
process X to another birth-and-death process X+ with birth rates b'^ = \n , ■ ■ ■ ,b~^ = Xi and 
zero death rates, in such a way that Xt/\TN — -^t' for all t > and X and X^ arrive in at 
the same time. In the present paper, we will extend their methods by showing that X and 
X^ can in addition be coupled to a process X~ with birth rates bf = Ai, . . . ,6^ = Xn and 
zero death rates, in such a way that X^ < Xiatjs, < ^'^^ all t > and all three processes 
arrive in N at the same time. 

1.2 Intertwining of Markov processes 

The coupling technique used by Diaconis and Miclo in |DM09] is of a special kind, which is 
sometimes called intertwining of Markov processes. Let X and X' be continuous-time Markov 
processes with finite state spaces S and S' and generators G and G', respectively, and let K 
be a probability kernel from S to S'. Then K defines a linear operator from M"^ to M.^ , also 
denoted by iT, by the formula 

Kf{x):=Y,K{x,y)f{y). (1.2) 

The following result, which is based on an observation by Rogers and Pitman |RP81) . was 
proved by Fill in |Fil92[ Thm. 2]. (An independent proof can be found in |AS101 Prop. 4]). 

Proposition 1.2 (Intertwining of Markov processes) Assume that 

GK = KG'. (1.3) 

Then there exists a generator G of an S x S' -valued Markov process with the property that if 
(X, X') evolves according to G and satisfies 

¥[X'^ = y\Xo\= K{Xo,y) {y G 5'), (1.4) 

then 

nx't = y I {Xs)^<s<t] = K{Xt, y) (t > 0, y G S'), (1.5) 

and the processes X and X' , on their own, are Markov processes evolving according to the 
generators G and G' , respectively. 

Algebraic relations of the type (jl.3p are called intertwining relations, hence the name inter- 
twining of Markov processes. We note that the operator K needs in general not have an 
inverse, and even if it does, this inverse will in general not be associated to a probability 
kernel from S' to S. In view of this, an intertwining of Markov processes is not a symmetric 
relation. To express this, following terminology introduced in [AS 10] . we will also say that in 
the set-up of Proposition 11.21 X' is an averaged Markov process on X. 

1.3 Intertwining of birth-and-death processes 

We are now ready to formulate our main result. Deviating slightly from our notation in Sec- 
tion [TTTl we let X = {Xt)t>Q be a continuous-time Markov process with state space {0, . . . , A^}, 
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started from Xq = 0, that jumps from x — 1 to x with birth rate bx and from x to x — 1 with 
death rate dx, where bi,. . . ,bN > 0, di,. . . ,d]\f-i > 0, but djy = 0, i.e., X is the stopped 
process from Section [TTTl We let G denote the generator of X, i.e., 

Gf{x):=bx+i{f{x + l)-f{x))+dx{f{x-l)-f{x)) {0<x<N), (1.6) 

where / : {0, . . . , N} — )• M is a real function and we adopt the convention that do = and 
bN+i = so that the corresponding terms in (II. 6p are zero, regardless of the (fictive) values 
of / in —1 and iV + 1. The following theorem is our main result. 

Theorem 1.3 (Intertwining of birth- and- death processes) The operator G has + 1 
distinct eigenvalues = — Aq > — Ai > • • • > —Xn- Let X^ and X^ he the pure birth processes 
in {0, ...,A'^}, started from Xq = Xq = 0, with birth rates b^ := Ai,...,6^ := Ajv and 
bf := Xn, . . . , 6^ := Ai, respectively, and let G^ and G'^ be their generators. Then there exist 
probability kernels K~ and on {0, . . . , N} satisfying 



K-{x,{0,...,x}) = l, K+{x,{0,...,x}) = l, iO<x<N) 
K-{N,N) = 1, K+{N,N) = 1, 

and 



[1.7) 



(i) K+G = G+K+ and (ii) GR- = R-G'. (1. 
Moreover, the processes X~,X, and X^ can be coupled in such a way that 

(i) nXt = y I (X+)o<.<i] =K+(X+,y) {t>0, 0<y<N), 

(ii) F[Xi = y 1 {Xt,Xs)o<s<t] = R-{Xt,y) {t>0, 0<y<N). 



(1.9) 



The existence of a kernel R'^ such that (jl.Sp (i) and (jl.9p (i) hold has been proved before 
in [DM09t Prop. 10]. Our new contribution is the construction of the kernel R~ such that 
moreover (jl.Sp (ii) and ()1.9p (ii) hold. It is easy to see that formulas ()1.7p and (|1.9p imply 
that 

(i) X7<Xt<X+ (t>0), 

, (1.10) 

(ii) T'j^ = TN = r]f, 

where r^v := inf{t > : Xt = N} and and are defined similarly for X~ and X^, 
respectively. We note that X~ and X~^ move, in a sense, in the slowest resp. fastest possible 
way from to N, given that they have to arrive at exactly the same time as X. Note that, 
using terminology introduced at the end of Section 11.21 ^ is an averaged Markov process on 
X~^ and X~ is an averaged Markov process on X. 



1.4 Discussion 

In comparison to the paper by Diaconis and Miclo [DM09], the present paper does not add too 
much that is new. In particular the construction of the kernel R~ in Theorem ll.3l is very similar 
to the construction of the kernel R~^, which was already carried out in [DM09] . However, we 
believe that the observation that both constructions are possible, with an interesting symmetry 
between them, is of some interest. 

The (new) construction with the process X~ has in fact one advantage over the construc- 
tion with X~^ , since Proposition [L2] and formula (II. Sp imply that the process X started in any 
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initial state can be coupled to a process X~ with the same dynamics as in Theorem 11.3^ in 
such a way that P[^(~ = y \ {Xs)Q<s<t\ = K~{Xt,y) for all < y < and t > 0. This implies 
that for a general initial state Xq = x £ {0, . . . ,N}, the stopping time Tjy is distributed as 
'^y=z where di , . . . , gn are independent exponentially distributed random variables with 
parameters Ai, . . . , Aat and Z is an independent {0, . . . , A^}-valued random variable with law 
K^{x, •). Note that the (old) coupling with the process X+ forces one to start the process 
X in an initial law that is a convex combination of the laws K~^(x, ■ ) with < j; < A'^, hence 
no conclusions can be drawn for arbitrary initial states. 

On the other hand, the methods of [DM09] can also be used to study birth-and-death 
processes on {0, . . . , N} whose death rate is not zero and which, therefore, converge in 
law to a unique equilibrium. In particular, Diaconis and Miclo use a generalization of their 
intertwining relation (jl.Sp (i) to construct a fastest strong stationary time for such processes 
(we refer to |DM09j for the definition). In contrast, it seems that the interwining relation 
(jl.Sp (ii) does not generalize to such a setting. 

On a more general level, one may ask what the advantage is of a 'probabilistic' proof of 
Proposition 11.11 as opposed to older, more algebraic proofs. Since most of the work behind 
Theorem 1 1 . 31 goes into proving the intertwining relations (jl.Sp . one might even argue that the 
present proof is still rather algebraic in nature, although with a strong probabilistic flavour. 
In this context, it is interesting to note that the fact that G is diagonalizable with real, 
distinct eigenvalues follows as a result of our proofs (in particular, this follows from a repeated 
application of the Perron-Frobenius theorem) and does not have to be provided by some extra 
argument (based on, for example, reversibility). 

In general, diagonalizing a generator of a Markov process gives very strong information 
about the process, but in practice, if the state space is large, it is hard to get good information 
about the position of eigenvalues etc. The idea of interwining generators with transition 
kernels may in some cases be a good way to transform generators of complicated processes 
into generators of more simple processes and thus provide a more probabilistic alternative to 
diagonalization. 

The methods of this paper can certainly be extended to one-dimensional processes with 
two traps, to dicrete-time processes, and to one-dimensional diffusions. Miclo |MiclO| has 
proved a generalization of Proposition fTD for reversible Markov chains. In [AS 10], intertwining 
relations were used to estimate the time to extinction for large hierarchical contact processes. 
The present work was partly motivated by an open problem from that paper. (To be precise. 
Question 1° from Section 3.3.) 

2 Proofs 

2.1 Leading eigenvectors 

Let X be the birth-and-death process in Sn := {0, ...,A^} from Section 11.31 and let G : 
jg'^AT _^ jg-Siv ]-,g j^g generator, defined in ()1.6p . We equip M'^'^ with the usual inner product 
(vrl/) := X^l^o ''^i^)fi^) denote the adjoint of G with respect to this inner product. 

Then 

G^tt{x) = bxTT{x - 1) - bx+iTr{x) + dx.+i7r(x + 1) - (ix7r(x), (2.1) 

where as in ()1.6p we use the convention that do = and &Ar+i = so that the corresponding 
terms in ()2.ip are zero, regardless of the (fictive) values of vr in — 1 and A^ + 1. 
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Since 5^ (the delta mass in A^) is the unique invariant law of X, the eigenvalue of 
the generator G has multiplicity one and its unique left and right eigenvectors are 5]sf and 
the constant function 1, respectively. We will need the following result on the next largest 
eigenvalue and its left and right eigenvectors. 

Lemma 2.1 (Leading eigenvectors) There exists a A > and fjir^ M'^'^ such that 

(i) / is strictly decreasing on {0, . . . , N} and satisfies /(O) = 1, f{N) = 0, 

(ii) vr is strictly positive on {0, . . . , — 1} and satisfies "^^=0 ^(^) = 1 = — 7r(A^), 

(iii) Gf = -Xf and G^ir = -Avr. 

Proof Set 

e{x) ■.= 5x {Q<x<N -I) 

{{x) ■=6^-6n lo<x <N - 1) 

Then {e(0), . . . , e{N)} is a basis for M^^ and {^(0), . . 
(e(x)l^(y)) = l{x=y}- Set 

£: := span{e(0), . . . , e(iV - 1)} = {/ G M^^ : /(iV) = 0}, 

:= span{e(0), . . . , ^(iV - 1)} = {vr G M^^ : ^(^) = 0}- ^'^ 

Since is a trap for the process X, it is easy to see that the operator G maps the space £ 
into itself. Since the coordinates of a vector in £ with respect to the basis {e(0), . . . ,e(A^)} 
are the same as its coordinates with respect to the standard basis {6o, . . . ,6]y}, it follows that 
with respect to the basis {e(0), . . . , e(A^)}, the matrix [G] of G has the form 




where A{x,y) = G{x,y) for < x,y < — 1. The restriction of the process X to the space 
{0, . . . ,N — 1} is irreducible in the sense that there is a positive probability of going from 
any state to any other state. Therefore, by applying the Perron-Probenius theorem (see, e.g. 
Chapter XIII, §2, Theorem 2 in |GanOO| ) to ^4 + cl and its adjoint for some sufficiently large 
c, one finds that A has a real eigenvalue —A of multiplicity one, which is larger than all other 
real eigenvalues, and associated left and right eigenvectors vr G J-" and f £ £ that are strictly 
positive with respect to the bases {■^(O), . . . , ^(A^ — 1)} and {e(0), . . . , e{N — 1)}, respectively. 
Since Markov semigroups are contractive we have — A < and since the eigenvalue zero of 
G has multiplicity one and belongs to different left and right eigenvectors, we conclude that 
—A < 0. Since we can always normalize our eigenvectors such that "^^^q Tr{x) = 1 and 
meiX^~Q f{x) = 1, this proves all statements of the lemma except for the fact that / is strictly 
decreasing. 

To prove this latter fact, we observe that by the facts that Gf = —Xf and / > on 
{0,...,N-1}, 

6i(/(l)-/(0))=-A/(0)<0, (2.5) 
which show that /(O) > /(I). By the same argument, 

6,+i(/(x + l)-/(x)) =-A/(x)-4(/(x-l)-/(x)) <0 {l<x<N-l), (2.6) 

from which we see by induction that f{x) > f{x + 1) for all < a; < A^ — 1. I 



and e(A^) := 1, , . 

and ^(N) := 6n. ^ ' 

, ^(A^)} is its associated dual basis, i.e.. 
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2.2 Intertwining the fast process 

In this section, we prove the existence of a kernel satisfying (ll.7p and (ll.Sp . Our proof 
is basically the same as the proof given in [DM09] . but as a preparation for the next section 
it will be convenient to review their proof and shorten it somewhat. The proof in |DM09] is 
written in such a way as to make clear how the authors arrived at their argument and uses 
discrete derivatives that are presumably also useful if one wants to generalize the theory to 
one-dimensional diffusions. If our only aim is Theorem 11.3] however, we can summarize their 
arguments quite a bit. 

The kernel will be constructed as the concatenation of an inductively defined sequence 
of kernels K^^^^^^ , . . . ,K^^^^ . Associated with these kernels is a sequence of generators 
. . . , + of birth-and-death processes in {0, . . . , A^} satisfying the intertwining 

relations 

^(M)+^(M)+ ^ ^(Af»l)+^(Af)+ (1 < M < iV - 1), (2.7) 

where the process with generator G^^^^ has birth rates 6^^^^ , • • • , b^lf "^ > and death rates 
d[^\ . . . ,d^^P > 0, d^^^^^ = . . . = (i^^^ = 0; see Figured] for a picture. In particular, we will 
choose G^^~^^~^ := G and setting G~^ := G^^^~^ will yield the desired pure birth process with 
birth rates 6^ = Aat, . . . , 6^ = Ai. 

The core the proof is the following proposition, which corresponds to the inductive step in 
the argument. 

Proposition 2.2 (Inductive step) Let I < M < N — I and let G be the generator of a birth- 
and-death process in {0, . . . , N} with birth rates 6i, . . . , 6jv > and death rates di, . . . , dM > 0, 
dM+i = • • • = div = 0. Then there exists a probability kernel K on {0, . . . , satisfying 

K{x,{<d,... ,x]) = 1 {{)<x<N) and K{x,x) = l {M + l<x<N), (2.8) 

and a generator G' of a birth-and-death process in {0, . . . , A^} with birth rates b'l, . . . ,b'j^ > 
and death rates d[, . . . , d'j^j_^ > 0, d\^ = ■ ■ ■ = d'j^ = 0, such that KG = G'K . 

Proof It follows from Lemma 12.11 applied to the process stopped at M + 1 that there exists 
a function p : {0, . . . , A^} ^ M such that p > on {0, . . . , M}, p = on {M + 1, . . . , N}, 
Ylx=oPi^) = 1, and 

G^p{x) = -Xp{x) + X6m+i{x) {0<x<N), (2.9) 

where 

A = bM+ip{M) > 0. (2.10) 

The law p is sometimes called a quasi- stationary law. Using p, we define the kernel K on 
{0,...,Ar}by 

K{x,y) := { ^ H[x) (2.11) 

l{j/=^l if M + 1 < X, 



where 



H{x):=Y,p{y) (0<x<M). (2.12) 
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G 



^(2) + 



^(0) + 




Figure 1: Intertwining of birth and death processes. In this picture = 4. All nonzero 
transition rates and probabilities have been indicated with arrows. 
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Since K is a. lower triangular matrix, it is invertible, so there exists a unique linear operator G' 
satisfying KG = G'K and G' is in fact given by G' = KGK-\ Since G'l = G'Kl = KGl = 
we see that 

G'{x,x) = -Y,G'{x,y). (2.13) 

In view of this, to prove our claim, it suffices to check that the off-diagonal entries of G' 
coincide with those of a birth-and-death process in {0, . . . , A^} with birth rates b'^, . . . ,b'j^ > 
and death rates d[, . . . , d'j^_^ > 0, d\,j = ■ ■ ■ = d'j^ = 0. 

To determine the off-diagonal entries of G' , we calculate, using (12. ip and (12. 9p . 



{KG){x,y) = G^K{x, 

p{y) 



-Al 



-V(y) + A5Af+i(y) 



dx+i^^T7^Uy) + if X < M, 



Hix) 



H{x) 



(2.14) 



if X = M, 
if X > M. 



In order to find C, we need to express these formulas, as functions of y, as linear combinations 
of the basis vectors {K{x^ ■ ))o<x<n- To that aim, we observe that 



5x = K{x,-) (M + l<x<iV), 
while for 1 < x < M, we have 

p{y) 



(2.15) 



^x{y) 



'1 



{y<x} 



1 



H{x) 

p{x) 
H{x) 

p{x) 



{y<x} 



p{y) H{x - 1) 



Hix) 



p{x) 



p{y) 



H(x - 1) 



(2.16) 



K{x,y) 



H{x - 1) 

p{x) 



K{x-l,y). 



Inserting this into (j2.14p . we find that 



N 



Y,G'{x,x')K{x',y) := {KG){x,y) 



x'=0 



-XK{0,y)-di 



-K{0,y) 



H{0) 

p(0) fHil) HiO) 
-Ai^(x,,)-4.i%tl)(^(-) 



+bx^ 



i7(x) V p{x 
p{x) (H{x + 1 



K{x,y)-^^^^—^K{x-l,y) 
p{x) 



H{x) \ p{x + 1) 
-XK{M,y) + XK{M + l,y) 

-bx+iK{x, y) + bx+iK{x + 1, y) 



K{x + l,y)-^^^K{x,y) 



p(x + 1) 



if X = 0, 



if < X < M, 

if X = M, 
if X > M. 



(2.17) 
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From this, we can read off tlie off-diagonal entries of G' . Indeed 

p{x)H{x + 1) 



61 



x+l 



G'{x,x + 1) = < 



+1 



d', = G'ix,x-l) 



A 

bx+i 

dx+i 




H{x)p{x + 1) 



if X < M, 
if X = M, 
if X > M, 



(2.18) 



p{x + l)H{x - 1) 
H{x)p{x) 



if < X < M, 
if X > M, 



and all other off-diagonal entries are zero. I 

Remark The proof of Proposition 12.21 is straightforward except for the clever choice of K in 
(j2.1ip - (|2.12p . For some motivation of this choice and the way the authors arrived at it we 
refer to |DM09 j. 

Using Proposition 12.21 we can construct a sequence of kernels K^^^'^^^ , . . . ,K^'^'>^ and gener- 
ators G*-^"^^'"*", . . . ,G^'^'>^ satisfying the intertwining relations (|2.7|) . such that G~^ := G^^^^ 
is a pure birth process with birth rates bf , . . . > 0, say. It is now easy to see that the 
composed kernel 

^+ .^^(i)+...^{JV-i)+ (2.19) 
satisfies K+{x, {0, . . . , x}) = 1 (0 < x < A^), K+(N, N) = 1 and K+G = G+ K+ . It is straight- 



forward to check that the eigenvalues of G' are —6 



1 ' ' 



-6^, 0. Since G 



{K+)-'^G+K+, 

the operators G and G~^ have the same spectrum. 

We claim that b^ > ■ ■ ■ > b^ > {). To see this, recall from the proofs of Lemma 12.11 and 
Proposition 12.21 that — fej^j is the Perron-Frobenius eigenvalue of the process with generator 
q{M)+ stopped at M -|- 1. It follows from the intertwining relation (12. 7p that —b\,j_^ is also 
an eigenvalue of this process, corresponding to a different eigenvector, hence by the Perron- 
Frobenius theorem, 6a/-i > 6m- 



2.3 Intertwining the slow process 

In the previous section, we have constructed a kernel and generator of a pure birth process 
G"*" such that (II. Sp (i) holds. In this section, we construct a kernel K~ and generator of a 
pure birth process G~ satisfying (jl.Sp (ii). The proof will be very similar to the previous 
case, except that some things will 'go he other way around'. In particular, using terminology 
introduced at the end of Section [1.21 G~ will be the generator of an avaraged Markov process 
X~ on X while in the previous section we constructed a pure birth process such that X 
is an averaged Markov process on X^ . 

As in the previous section, the kernel K~ will be constructed as the concatention of an 
inductively defined sequence of kernels i^*^^^ ~, . . . , iT^^"^-* ~. Associated with these kernels 
is a sequence of generators G^^-*'", . . . , G^^"^) ~ of birth-and-death processes in {0, ...,A/^} 
satisfying the intertwining relations 

^(M-l)-^(M)- ^^(M)-^(M)- (1<M<7V-1), (2.20) 

where the process with generator G^^^^ has birth rates 6^^^^ , • • • , 6;^^^ > and death rates 
^(M) _ _ , _ ^{M) _ d^^^^, . . . > 0, and d^^^^ = 0. We again refer to Figured] for an 

illustration. 

The core of the argument is the following proposition. 
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Proposition 2.3 (Inductive step) Let < M < N — 2 and let G he the generator of a 
hirth- and- death process in {0, . . . , A^} with birth rates 6i,...,6jv > and death rates di = 
■ ■ ■ = dM = 0, djvf+i; • • • ) djv-i > 0, and dj^ = 0. Then there exists a probability kernel K on 
{0, . . . , A^} satisfying 

ir(x,{0, ... = 1 (0<x<iV) and K{x,x) = l (x {M, . . . , iV - 1}), (2.21) 

and a generator G of a hirth- and- death process in {0, . . . , N} with birth rates 6i, . . . , 6jv > 
and death rates di = ■ ■ ■ = dM+i = 0, dM+2, ■ ■ ■ , d^-i > 0, and d^ = 0, such that GK = KG. 

Proof It follows from Lemma |2 . 1 1 applied to the process restricted to {M, . . . , A^} that there 
exists a function / : {0, . . . , N} M such that / = on {0, . . . , M — 1}, / is strictly decreasing 
on {M, . . . , TV}, /(M) = 1, f{N) = 0, and 



where 
We set 



Gf{x) = -\f{x) + hMh.i-i{x) (0 < X < iV), 
A = 6m+i(1-/(M + 1)) >0. 



(2.22) 
(2.23) 



K{x,y):=l{^=y} {y^{M,...,N -I}). (2.24) 

For y = M, . . . ,N — 1, we claim that we can inductively define the kernel K{x, y) and contants 
Gy > in such a way that 



(i) K{x,y) := Cyl{y<^yf{x), 

(ii) i^(y,y') = l, 

y'=M 



{M <y<N -I) 



(2.25) 



To see that this is all right, note that for y = M (j2.25p (i) and (ii) are satisfied by choosing 
Gm '■= 1; while for Af + 1 <y<A^ — 1 (I2.25P (i) and (ii) imply that we must choose 



p.-7^(i-E'^-fe.!''))- 

■' y'=M 

Since / is strictly decreasing on {M, . . . , A^}, one has, by induction. 



(2.26) 



^ K{y,y')= Cy>f{y)< Gy>f{y-l)= K{y-l,y') = l, 



(2.27) 



y'=M 



y'=M 



y'=M 



y'=M 



which shows that Gy > 0. We now calculate 
{GK){x,y)=GK{-,y)ix) 

hy5y-l{x) - by^l5y{x) 

-\f{x) + hMhi-i{x) 

-\Gyl{y<.j,}f{x) + hyGy5y-.l{x) 

-dyGy{f{y-l)-f{y))5y{x) 



if < y < M - 1, 
if y = M, 

ifM + l<y<A^-l, 
iiy = N. 



(2.28) 
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By the same arguments as those in the previous section, there exists a unique hnear operator 
G such that GK = KG. In order to check that G is the generator of a birth-and-death 
process in {0, ... , A^} with birth rates 61, . . . , 6jv > and death rates di = ■ ■ ■ = dM+i = 0, 
(iM+2) • • • ) d^-i > 0, and djy = 0, it suffices to check that the off-diagonal entries G{x, y) have 
the desired form. In order to do this, we must express the formulas in ()2.28p . as functions of 
X, as hnear combinations of the basis vectors {K( ■ ,y))o<y<N ■ We observe that 



6y{x) = K{-,y) 
while for M < y < N — 2, we have 

6y{x) = {l{y<x} - '^{y+l<x}) 



1 



(y0{M,...,iV-l}), 

M 
m 



(2.29) 



f{y)Cy 



Cyl{y<x}fix) 



1 



and 



f{y)Cy 



K{x,y) 



1 



f{y)Cy+i 



Cy+l^{y+l<x}f{x) 



(2.30) 



f{y)Cy+i 

1 



f{N - 1)Gn- 



K{x,y + l), 
-K{x,N - 1). 



(2.311 



Inserting this into (I2.28p . we obtain 

Y,Kix,y')G{y',y) = {GK)ix,y) 
y' 

byK{x, y -1) - by+iK{x, y) 
-XK{x, M) + hMK{x, M-l) 



K{x,y) 



-Ai^(-,.) + ^.C^.(^^_^ /(, + l)a 
-d,C,(/(,-l)-/(,))^^(^'^) ^(^'^ + 1 



f{y)Gy /(y)C,+i 
-dyCy{f{y-l)-f{y)f^''^y^ 



if < y < M - 1, 
if y = M, 



if M + 1 < y < iV - 2, 



6iv 



K{x,N -I) 
f{N - 1)Cn^i 



mcy 



ify = N-l, 
if y = AT, 



(2.32) 

where we use the convention that bo = and hence boK^x, —1) =0, regardless of the (fictive) 
value of K{x, —1). From ()2.32p we can read off the off-diagonal entries of G. Indeed, 



G(y-l,y) 



4+1 = G{y + i,y) ^ 



'f{y-l)Cy-i 



if 1 < y < M, 

if M 1 < y < A^, 



C,(/(y-l)-/(y)) 
Cy+if{y) 



if y^{M + 1,..., AT -2}, 
if M -Fl < y < Af - 2, 



(2.33) 
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and all other off-diagonal entries are zero. We note that in particular, by (j2.23p and the 
definition of the C^'s, 

bM+i = ^yA'n^^' = ^m+iCm+i = bM+i (1 - /(M + 1)) = A. (2.34) 



Remark As in the case of Proposition 12.21 the proof of Proposition 12.31 is straightforward 
except for the choice of the kernel K. We have guessed formula ()2.25p by analogy with 
formula (f2ll]) . which is due to |DM09j . 

With the help of Proposition 12.31 we can inductively define kernels K^^^^, . . . ,K^^^^^^ and 
operators G^^)'", . . . , Setting G" := G^^-^) " and 

i^- = i^«----K(^-i)- (2.35) 

now yields a generator of a pure birth process with birth rates b'^ , . . . and a kernel K~ 
with the properties described in (jl.7p - (|1.8p . 

In the same way as in the previous section, we see that 0, — 6]", . . . , —bjf are the eigenvalues 
of G. To see that < 6i < • • • < ^at we observe from (j2.34p that — fej^f Perron-Probenius 
eigenvalue of the process with generator 

q(M) - i-estricted to {M, ...,N}. It follows from the 
intertwining relation (j2.20p that — is also an eigenvalue of this process, corresponding to 
a different eigenvector, hence by the Perron- Frobenius theorem, bM < &a/+i- 

2.4 Proof of the main theorem 

Proof of Theorem 11.31 The existence of generators G^,G^ and kernels K~,K^ satisfying 
(jl.7p - (|1.8p has been proved in the previous sections. By Proposition 11.21 it follows that 
and X can be coupled such that (jl.Op (i) holds. By applying Proposition 11.21 to the kernel L 
from {0, . . . , N}'^ to {0, ... , A^} given by 

L((x, y), z) := ^-(y, z) (0 < x, y, z < N), (2.36) 

we see that (A+, A) and can be coupled in such a way that both (jl.9p (i) and (ii) hold.I 
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